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Abstract. Asteroseismology is found to be a excellent tool for detecting diffusion- 
induced helium gradients inside main-sequence A stars. Models have been com- 
puted for 1.6 and 2.0 Mq stars with pure helium diffusion, at different ages, so 
that the helium gradient lies at different depths inside the star. The adiabatic 
oscillation frequencies have been analysed and compared with those of a model 
without diffusion. Clear signatures of the diffusion-induced helium gradient are 
found in the so-called "second differences" : these frequency differences present 
modulations due to the partial reflexion of the sound waves on the layer where 
the helium gradient occurs. A tentative application to the roAp star HD60435, 
which presents enough detected oscillation frequencies for the test to be possible, 
is very encouraging. The results suggest the presence of a helium gradient inside 
the star, which is consistent with the idea that the triggering of the oscillations is 
due to the hydrogen K-mechanism. 
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1. Introduction 

It is now widely recognized that the abundance anomalies observed in peculiar A stars 
are basically due to element diffusion. Those for which the radiative acceleration is larger 
than gravity (like most of the metals) are pushed upwards while the others, such as helium, 
diffuse downwards fMichaud ll97UI Vauclair & Vauclair II 982(1 . The resulting atmospheric 
Send offprint requests to: S. Vauclair 
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abundances depend on the competition between element diffusion and macroscopic mo- 
tions like convection, mass loss, rotation-induced mixing and other kinds of mixing pro- 
cesses. Among the peculiar A stars, those which have large organized magnetic fields 
(the so-called Ap stars) are more complex than the others (like Hg-Mn stars) : they show 
evidence of abundance spots or rings modulated by the magnetic structure. Among the 
coolest Ap stars (effective temperatures between 6700K and 8700K) some oscillate with 
periods of a few minutes (the roAp stars) while the others seem stable (noAp stars). In 
roAp stars, the amplitudes of the oscillation modes are modulated according to the ro- 
tation period, which is explained in the framework of the oblique pulsator model (Kurtz 

The helium K-mechanism was invoked in the past as a possible way to trigger the 
oscillations : a well-adjusted competition between helium settling and mass loss could 
increase the helium abundance in stellar atmospheres as observed in helium-rich stars 
fVauclair I1975L Vauclair, Dolez & Gough [l991(l . However, recent papers (Dziembowski 
& Goode [T996l Balmforth et al. l2[)t)l|l show that, even in this case, helium is not able to 
destabilize the star. On the other hand they find that, in some cases, the oscillations in 
A stars can be triggered by the hydrogen K-mechanism. Such a process may be enhanced 
by the settling of helium, which diffuses downwards in the absence of strong mass loss. 

In the present paper we study the asteroseismic signature of a helium gradient inside 
main-sequence A- type stars. As mentioned by Gough l)1990|l . rapid variations of the sound 
velocity inside a star lead to partial reflections of the sound waves, which may clearly 
appear as frequency modulations in the so-called "second differences" : 

^2t'n,i = i^(„+l),Z + i^(ri-l),; ~ 2i/„^; (1) 

which are computed for the same value of the azimutal number I. The modulation period 
of the oscillations is twice the "acoustic depth" of the region where the feature occurs 
(i.e. the time needed for the sound waves to travel between this region and the stellar 
surface). 

Such an effect has been extensively studied in the literature for stellar models in which 
helium settling is neglected. In this case, the modulations in the frequencies are mostly 
due to the Hell ionisation zone and to the edge between the convective and radiative 
zones (Monteiro & Thompson 119981 Roxburgh & Vorontsov 120011 Mazumdar & Antia 
IWn Migho et al.EnSni) 

We show here that the presence of a diffusion-induced helium gradient leads to a 
kink in the sound velocity with a very clear signature in the oscillation frequencies : 
asteroseismic observations of A-type stars can test helium diffusion and lead to a precise 
value of the acoustic depth corresponding to the position of the helium gradient inside 
the star. 
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For the present study, we neglect the effects of magnetic fields (the stars are assumed 
spherically symetric) and we do not compute the diffusion of heavy elements. Complete 
computation of element diffusion including the radiative acceleration on metals is a de- 
manding work, especially for these stars where heavy element diffusion occurs in the 
atmosphere, so that radiative transfer should be solved in detail in the optically thin 
regions. This is out of the scope of the present paper, where we are interested in the 
structural changes induced by helium settling below the stellar surface and its influence 
on the oscillation frequencies. It has recently been shown (e.g. Richard, Michaud & Richer 
I2001|l that in case of pure diffusion iron should accumulate deeper inside the stars, at 
temperatures around 200 000 K, where it can create a small convective zone. As the 
results we obtain here are very encouraging, new computations including heavy element 
diffusion should be done in the near future to see whether such a feature may be checked 
with asteroseismology. 

We apply this test to the only roAp star in which enough modes have been observed, 
namely HD 60435 (Matthews et al. I1987|l . The results are extremely encouraging and 
show the importance of detecting as many modes as possible in rapidly oscillating Ap 
stars. 

2. Theoretical discussion 

2.1. Seismic signatures of helium gradients 

Stellar acoustic p modes with low I values can propagate deeply inside the stars. For this 
reason, they may be used to obtain information on the deep stellar structure. However, in 
the case of strong gradients in the sound velocity, which may be due to the boundary of 
a convective zone, to the helium ionization region, or to helium gradients, the waves are 
partly reflected : this creates modulations in the frequency values which may be clearly 
visible in the second differences. The modulation periods are equal to 2ts where ts is 
the time needed for the acoustic waves to travel between the surface and the considered 
region (acoustic depth), i.e. : 



where c(r) is the sound velocity at radius r, and the radius of the considered region. 

As the I dependence of the oscillatory signal is small for high frequency modes of low 
degrees (Monteiro, Christensen-Dalsgaard & Thompson ll994p . computational results for 
modes of different low I values can be treated together. Here we show that the presence 
of diffusion-induced helium gradients lead to very clear signatures in the modulation of 
the second differences. 




(2) 
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TABLE 1 
Model parameters 



mass (M0) 


age (Myrs) 


log(L/L0) 




tx (sec) 


ts (sec) 


Ts/R 


1.6 


95 


0.84 


7936 


5400 


1480 


0.93 


1.6 


591 


0.88 


7585 


6920 


2700 


0.86 


1.6 


1600 


1.05 


7071 


11240 


5160 


0.80 


1.6h 


1600 


1.06 


7352 


10227 






2.0 


63 


1.22 


9521 


5917 


1470 


0.94 


2.0 


649 


1.31 


8293 


10150 


4100 


0.85 



Note.- tx represents the total acoustic radius of the models, ts the acoustic depth at the location 
of the helium gradients and Vg/R the corresponding fractional radii ; all the models presented 
in this table include helium diffusion, except the model labelled 1.6h which is homogeneous in 
its outer layers. 

2.2. Models with helium diffusion 

We have chosen to study the evolution of 1.6 M0 and 2.0 M© stellar models, with 
pure helium settling. The stellar evolution code we use is the "Toulouse- Geneva code" 
which has been described many times in the literature (e.g. Richard et al. 119961 Theado 
& Vauclair I2003|l . The physical input and parameters are those described in Richard, 
Theado & Vauclair l|2()()4|l : they include the most recent studies available for the equation 
of state, opacities and nuclear reaction rates. Helium diffusion is treated as described in 
these papers (see also Vauclair l2U(J3f) . Metal diffusion is not included here and no mixing 
process or mass loss is taken into account. In real stars, there must be some mixing which 
prevents element diffusion from producing extreme abundance anomalies. Computations 
including mixing as well as metal diffusion will be done in the near future. The aim of the 
present paper is to study the precise signature on the oscillation frequencies of helium 
gradients inside stars. 

We have computed the oscillation frequencies and their second differences in five 
models : 1.6 Mq with ages 95 Myrs and 591 Myrs, and 1.6 Gyrs and 2.0 M© with ages 63 
Myrs and 649 Myrs (Table 1). The results are given in Figs. 1 to 5. Only modes of degrees 
/ = 0, 1, 2 and 3 have been taken into account so that they could be treated together for 
the study of the second differences. The oscillation frequencies have been limited to the 
range 0.5 to 2 mHz to take into account the asymptotic approximation validity and the 
stellar cut-off. In this range 80 to 100 second differences are computed, according to the 
model. 

In each figure six graphs are displayed : a) the helium profile inside the star ; b) the 
sound velocity, which shows a clear kink at the place of the helium gradient ; c) the first 
derivative of the sound velocity in which the kink is still clearer ; d) the second differences 
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which show periodical osciUations ; e) the Fourier transform of these osciUations, in which 
clear peaks are found for precise time values ; f) the time needed for the acoustic waves 
to travel between the surface and the considered radius, or "acoustic depth" . 

We obtain clear signatures in the second differences of the frequencies of the diffusion- 
induced helium gradients inside stars. It can be checked that the periods corresponding 
to the peaks in the Fourier transforms (graphs e) are exactly twice the acoustic depths 
of the helium gradients (graphs f), as expected. In these figures, the scales are all the 
same except for graphs d) and f). In particular the Fourier transforms (graphs e) are 
presented at the same scale in each figure to show how the peak amplitudes decrease 
for deeper layers (larger times). As shown by Mazumdar & Antia (|2001(l . the amplitude 
of the oscillatory signal in the second differences contains an amplification factor of 
4sin^(7rti/t,) where t* is the total acoustic radius of the star and ts the acoustic depth 
of the partial reflection region. Here this factor increases from about 2 to about 4 as 
the models evolve. The amplitude decrease by a factor ~ 10 obtained in the present 
computations must be due to a competing process. It may be due partly to the fact that 
the energy in p waves decreases towards the center and partly to the fact that the helium 
gradient becomes smoother with age, so that the reflected part (oscillations) is smaller 
for deeper discontinuities. 

With these figures we can clearly follow the signatures of the helium gradients as they 
sink into the stars. The corresponding ages are obtained with the assumption of pure 
diffusion. In real stars macroscopic effects may change the time scales for helium settling. 
Comparisons with observed frequencies will give the positions of helium gradients inside 
the star, but we must keep in mind that the time needed for these gradients to develop 
depends on the physics. 

2.3. Models without element diffusion 

Similar computations have been done for evolution models in which helium diffusion is 
suppressed. In this case no helium gradient exists in the stellar outer layers and the 
oscillatory signal described in the previous section is not present. However, helium ion- 
ization induces a typical variation of the Fi coefficient which gives rise to another type 
of oscillatory signal in the frequencies. Computational results are shown in Fig. 6 for the 
case of the 1.6 Mq model at 1.6 Gyrs. In this case graph a) displays the Fi coefficient in 
the stellar outer layers. As already studied by several authors (e.g. Miglio et al. I2003|l . 
this coefficient is strongly influenced by the helium ionisation gradients. We can see the 
corresponding kinks in the sound velocity and its first derivative. This oscillatory signal 
does not appear when diffusion is effective as helium is then depleted in the region con- 
sidered. Note that in the graphs of Fig. 6 the scales are quite different from those of Figs. 
1 to 5. The effect due to the helium ionization gradients is much smaller than the effect 
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Fig. 1. Helium diffusion and its consequences for the stellar structure and oscillation 
frequencies for a 1.6 Mq at 95 Myrs ; a) helium profile as a function of the fractional 
radius ; b) sound velocity : a kink is clearly visible at the place of the helium gradient 
; c) first derivative of the sound velocity : the kink is still more visible ; d) the second 
differences of the oscillation frequencies plotted as a function of the frequencies; e) the 
Fourier transform of graph (d) plotted as a function of time (in thousands of seconds) 
; the vertical dashed line corresponds to twice the total acoustic depth of the star (see 
Table 1) ; f ) the "acoustic depth" , or time needed for the acoustic waves to travel from the 
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Fig. 2. Same as Fig. 1, for a 1.6 M© model at 591 Myrs. 



due to the overall helium gradient in the case of diffusion. Here the amplification factor 
4sin^(7rfs/i*) is about 10 times smaller for the homogeneous model than for the model 
of the same age which includes helium diffusion. The large difference between the signal 
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Fig. 3. Same as Fig. 1, for a 1.6 Mq model at 1.6 Gyrs. 



amplitudes in the two cases may be due to this effect. The right-hand graphs display the 
second differences, their Fourier transform and the travel time of the waves, as in Figs. 1 
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Fig. 4. Same as Fig. 1, for a 2.0 Mq model at 63 Myrs. 



to 5. Tlie features are very different from those obtained in case of diffusion, where the 
signature of the helium gradient is dominant. 
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Fig. 5. Same as Fig. 1, for a 2.0 M© model at 649 Myrs. 
3. The roAp star HD60435 

Among all the rapidly oscillating Ap stars which have been detected up to now, only 
one, HD60435, presents a frequency spectrum rich enough for the second differences to 
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Fig. 6. This figure is similar to Fig. 1 wiien element diffusion is suppressed, in a 1.6 Mq 
at 95 Myrs ; here graph a) represents the Fi coefficient which shows a clear feature at 
the place of the Hell ionisation region ; the other graphs are similar to those of Fig. 1 ; 
we can clearly see the influence of the helium ionisation regions on the frequencies. Note 
that the effect is much smaller than that of a helium gradient when diffusion is taken 
into account. 
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be computed. This star has been extensively by Matthews et al (|1987|l who discovered 
17 frequency peaks in the fourier spectrum, which they basically identified as / = 1 and 
/ — 2 modes. As for all roAp stars, the amplitudes of the oscillations are modulated 
according to the rotation period (here 7.7 days). However, contrary to most Ap stars, no 
magnetic field has yet been clearly detected. 

We have computed the observed second differences for this star, using the frequencies 
given by these authors in their Table 2 (with an uncertainty of iO.OOOlmHz) and their 
tentative identification, as given in the same table. For the 0.8428 niHz frequency, we 
chose the first identification given, i.e (n — 12,2) and we assumed that the 1.2848 mHz 
frequency, which is not clearly identified in their paper, corresponds to (n — 3, 1). 

We present these frequencies and mode identifications in Table 2, with separate 
columns for I ~ 1 and 1 = 2. The second differences can only be computed when at 
least three consecutive modes have been identified for a given I. Table 3 gives the com- 
puted second differences for HD60435 and the I value of the corresponding mode. The 
obtained values for both I are presented in Fig. 7. The modes are labelled with their I 
value. 

The second differences clearly show a modulated trend similar to those presented in 
Figs 1 and 4, characteristic of the presence of a helium gradient below the surface. The 
fourier transform of the observed second differences is presented in Fig 8 (thick curve) 
together with those computed for 1.6 M© (95 Myrs) and 2.0 M© (63 Myrs). We can see 
that the observed peak for HD6G435 lies at a position close to those of the two models. 
The maximum is found at a period of 3400 sec, which corresponds to an acoustic depth 
of the helium gradient tg ~ 1700 sec. With a spectral classification A3V, this star should 
have a mass close to 2.0 Mq.Wc can deduce from this study that the modulation observed 
in the second differences for the star HD60435 may be a signature of a helium gradient 
at a radial depth of about 90% to 95% of the stellar radius. 

The uncertainties in the frequency measurements lead to relative uncertainties of 
the second differences smaller than 10%. Simulations of the error in the position of the 
amplitude peak due to the uncertainties in the frequencies show that the results are 
quite robust : the time tg at the maximum does not change by more than 0.05%. A 
more important uncertainty arises from the small number of observed frequencies. A 
reasonable estimate of this uncertainty is obtained by assuming that a new maximum 
could appear between the presently observed one and the two adjacent points if more 
frequencies were observed. Such an estimate leads to ±300 sec in the Fourier transform, 
oi ■.ts = 1700 ± 150 sec. 

These results relie on the mode identification which is given as "tentative" by 
Matthews et al ()1987|l . Note however that the difficulty in the identification was due 
to the fact that the frequencies were not exactly at the expected values. The present 
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v (mHz) mode id. : v (mHz) mode id. 



0.7090 


(n 


- 14, 1) 


0.8428 


(n 


- 12, 2) 


0.7614 


(n 


- 13, 1) 


0.9397 


(n 


- 10, 2) 








0.9906 


(n 


- 9, 2) 


1.1734 


(n 


- 5, 1) 


1.0433 


(n 


- 8, 2) 


1.2250 


(n 


-4, 1) 


1.0990 


(n 


- 7, 2) 


1.2848 


(n 


- 3, 1) 


1.1482 


(n 


- 6, 2) 


1.3281 


(n 


-2, 1) 


1.3525 


(n 


- 2, 2) 


1.3810 


(n 


- 1, 1) 


1.4073 


(n 


- 1, 2) 


1.4334 


(n 


, 1) 


1.4572 


(n 


, 2) 



Note.- after Matthews et al 1987 ; left columns -.1 = 1 right columns -.1 = 2. 

TABLE 3 

Computed second differences for HD60435 



V (mHz) sec.dif. (mHz) I 



0.9906 


+ 0.0018 


2 


1.0433 


+ 0.0030 


2 


1.0990 


- 0.0065 


2 


1.2250 


+ 0.0082 


1 


1.2848 


- 0.0016 


1 


1.3281 


+ 0.0096 


1 


1.3810 


- 0.0005 


1 


1.4073 


- 0.0049 


2 



study shows that the observed frequencies are consistent with the presence of a hehum 
gradient. 

4. Conclusion 

We have computed models of 1.6 and 2.0 main sequence stars with and without 
hehum diffusion and we have shown that asteroseismology can give clear signatures of 
the presence of helium gradients below the surface in these A-type stars. The analyses 
of the second differences in the oscillation frequencies show modulations which are due 
to the partial reflection of the acoustic waves at the place of the helium gradient. 
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At the present time these theoretical results can be compared to the observations of 
only one star, namely the roAp HD60435. The computations of the second differences 
can only been done if enough modes with adjacent values of n are detected for the same 
value of I. For HD60435, with the presently detected modes, we could compute 8 second 
differences, which enabled us to give evidence of a modulation similar to the one expected 
in case of the presence of a helium gradient below the surface. Furthermore, we could 
derive the acoustic depth of this gradient : 1700 ± 150 sec. The corresponding fractional 
radius is model-dependent, but from our computations we suggest that it should lie 
around r/R ~ 0.9 to 0.95. 

In our evolutionary sequences, which are computed with pure microscopic diffusion, 
a helium gradient at this depth is obtained early in the stellar life (65 Gyrs). However, 
in real stars microscopic diffusion must compete with macroscopic motions so that the 
helium gradient can develop later in the evolution : our strong result for HD60435 consists 
in determining the acoustic depth of the helium gradient, not the age of the star. 

As for all roAp stars, the amplitude of the oscillations observed in this star is modu- 
lated according to the rotation period, of 7.7 days. This modulation is probably related 
to the presence of a magnetic field, although no clear detection have been obtained. 
Our theoretical computations have been done for spherically symetrical stars, without 
rotation-induced mixing and ignoring any magnetic effect. In the future, the influence of 
the magnetic structure of the stars will have to be studied in detail. New models will also 
be computed, including the abundance variations of heavy elements (with computations 
of radiative accelerations), which have been neglected here. 

The results we have obtained are consistent with the idea that helium diffusion does 
occur in roAp stars. It has been suggested in the past that the oscillations observed 
in these stars could be triggered by helium K-mechanism ; such a process could occur 
only if helium accumulated due to the combined effect of diffusion and a hypothetical 
stellar wind (Vauclair, Dolcz, Gough [T991|l . More recently Balmforth et al l|2()01|l showed 
that such a helium accumulation was unable to destabilise the stars, while the triggering 
could be due to the hydrogen k- mechanism, enhanced in the polar regions in case of 
pure helium diffusion. The present results are strongly in favor of this model. 

In the future, new models will be computed, which will include the abundance varia- 
tions of heavy elements (with computations of radiative accelerations) . The influence of 
the magnetic structure of the stars will also have to be studied in details. 

It will be interesting to try to detect as many modes as possible in roAp stars to make 
similar comparisons with other stars. Few modes are generally detected in these stars, 
probably due to mode trapping. New instruments and techniques may help in detecting 
other modes. The frequencies should be determined with a precision of order 0.1 fiHz for 
the second differences to be computed with a relative uncertainty of order 10% and the 
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method described here to be apphcable. HD60435 should also be studied in more detail 
with the most recent techniques to obtain a higher precision of the results. 
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Fig. 7. Second differences in the oscillation frequencies of HD60435 as found from 
Matthews et al. (1987) ; the points are labelled with their I values. 
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Fig. 8. Fourier transform of the second differences of HD60435 (Fig. 7) in heavy sohd 
line ; the dashed curve represents the computed fourier transforms of the 1.6 M©, 95 Myr 
model and the light solid curve represents that of the 2.0 M©, 63 Myr model. 



